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Abstract 

We compute the 0{al) inclusive cross section for the process pp tih in the Standard Model, at ^/sff = 14 TeV. 
The next-to-leading order corrections drastically reduce the renormalization and factorization scale dependence of the 
Born cross section and increase the total cross section for renormalization and factorization scales larger than mj. 
These corrections have important implications for models of new physics involving the top quark. 
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1. One of the major goals of the LHC is uncover¬ 
ing the mechanisms of electroweak symmetry breaking 
and the generation of fermion masses. In the Standard 
Model of particle physics, the masses of gauge bosons 
and fermions are generated by a single scalar field. Af¬ 
ter spontaneous symmetry breaking, a neutral CP-even 
Higgs boson, h, remains as a physical particle. The 
fermion masses then arise through couplings to the Higgs 
boson. In the Standard Model, this coupling is directly 
proportional to the fermion mass. Since the top quark is 
the most massive quark, its coupling to the Higgs boson 
is particularly sensitive to the underlying physics. 


The associated production of a Higgs boson with a tt 
pair at the LHC, pp —> tth, will play a very important 
role in the 115 GeV< Mh <140 GeV Higgs boson mass 
range, both for discovery and for precision measurements 
of the Higgs boson couplings. This process will provide 
a direct measurement of the top-quark Yukawa coupling 
and will be instrumental in determining ratios of Higgs 
boson couplings in a model independent way 01- Such 
measurements could help to distinguish a SM Higgs bo¬ 
son from more complex Higgs sectors, e.g., as predicted 
by supersymmetry, and shed light on the details of the 
generation of fermion masses. 


2. The inclusive total cross section for pp tth at 
0{al) can be written as: 


,{pp tth) = ^ 


1 


1 -f Si 


dxidx2 


( 1 ) 


■[Tf{xi,p)T^^{x2,R)d’^i^,^{xi,X2,li) + (1 ^ 2 )] 


where Tf is the NLO parton distribution function for 
parton z in a proton, defined at a factorization scale 
pLf=pL, and ct’nlo is the 0(a\) parton level total cross sec¬ 
tion for incoming partons i and j, made of the channels 
qq, gg tth, and {q,q)g tih{q,q), and renormalized 
at the scale pr which we also take to be pr = R- At the 
LHC, the dominant contribution is from the gluon-gluon 
initial state, although the other contributions cannot be 
neglected and are included in this calculation. 

The NLO parton-level total cross section, con¬ 

sists of the 0(a1) Born cross section, (j^o, and the 0{as) 
corrections to the Born cross section, Scd^Lo, including 
the effects of mass factorization. Scr’MLo contains virtual 
and real corrections to the parton-level tih production 
processes, qq — > tih and gg — > tih, and the tree-level 
{q,q)g initiated processes, (q,q)g tih{q,q), which are 
of the same order in as ■ It can be written as the sum of 
two terms: 


In order to interpret the evidence for tih production 
and the measurement of the tih coupling as a verifica¬ 
tion of the Standard Model or as a signal for new physics, 
it is necessary to have a precise prediction for the cross 
section. QGD corrections are expected to be important 
and are crucial in order to reduce the dependence of the 
cross section on the arbitrary renormalization and fac¬ 
torization scales. In this letter, we present the next-to- 
leading-order (NLO) QGD corrections to the total cross 
section for pp tih at the LHG. Our results are in very 
good agreement with those of Ref. Q within the statisti¬ 
cal errors. Results for the Fermilab Tevatron have been 
presented elsewhere [^ ^ ||. 


= K^irt + Kial = (2) 

d{PS^)M{ij ^ tih) + J d{PS/i)M{ij ^ tih-\-k) , 

where M{ij tih) and M{ij tih -\- k) (for k = g,q, q) 
are respectively the matrix elements squared for the 
0{al) 2^3 and 2^4 scattering processes averaged 
over the initial degrees of freedom and summed over the 
final ones, while ^(PS's) and d{PS 4 ) denote the inte¬ 
gration over the corresponding three/four particle phase 
space. The {q, q)g initial state contributes only to 

The main challenges in the calculation come from the 
presence in the virtual corrections of pentagon diagrams 
with several massive external and internal particles, and 
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from the computation of the real part in the presence of 
infrared singularities. 

3. The 0(as) virtual corrections to the tree level 
ij —> tth {ij = qq, gg) processes consist of self-energy, 
vertex, box, and pentagon diagrams. The calculation of 
the virtual corrections to the qq initial state is described 
in Ref. [|j. The basic method is to reduce each diagram 
to a sum of scalar integrals of the form, 

/ pTRd n TTirr—^2-2 t ; po = 0, (3) 

that may contain both ultraviolet (UV) and infrared (IR) 
singularities. The finite scalar integrals are evaluated by 
using the method described in Ref. and cross checked 
with the numerical package FF 0. The scalar integrals 
that exhibit UV and/or IR divergences are calculated 
analytically. Both the UV and IR divergences are ex¬ 
tracted by using dimensional regularization in d = 4 — 2e 
dimensions. The UV divergences are then removed by 
introducing a suitable set of counterterms, as described 
in detail in Ref. Q . The IR divergences are cancelled by 
the analogous singularities in the soft and collinear part 
of the real gluon emission cross section. 

The most difficult integrals arise from the IR diver¬ 
gent pentagon diagrams with several massive particles. 
In Refs. ill we calculated the pentagon scalar inte¬ 
grals as linear combinations of scalar box integrals using 
the method of Ref. |,| 1 - For the gg initiated process 
we also used the method of Ref. and found perfect 
agreement between the results of the two methods. The 
virtual corrections to the gg initiated process have an 
additional complication with respect to the qq case be¬ 
cause of the presence of pentagon tensor integrals with 
rank higher than one. Pentagon tensor integrals can give 
rise to numerical instabilities due to the dependence on 
inverse powers of the Gram determinant (GD): 


GD(pi +p2-^'^pi) = - 


i—3 


[s - ( 2 mt + Mh)'^] 
64 


i^h 4“ ~ 545)^ ~ -|- S45)]sS45 X 

sin^ 6*45 sin^ ^45 sin^ 0, (4) 


where the 3-particle phase space has been expressed in 
terms of a time-like invariant S45 = (p4 -l-ps)^, polar an¬ 
gles, 045,0 and azimuthal angles, (j) 45 ,(f>, and s = xiX 2 Sh 
is the partonic center-of-mass energy squared. As can 
be seen in Eq. (Q), the Gram determinant vanishes when 
two momenta become degenerate, i.e. at the boundaries 
of phase space. These are spurious divergences, which 
cause serious numerical difficulties. We use two meth¬ 
ods to overcome this problem and find mutual agreement 
within the statistical uncertainty of the phase space in¬ 
tegration: 

• Impose kinematic cuts to avoid the phase space re¬ 
gions where the Gram determinant vanishes. Then 


apply an extrapolation procedure from the numer¬ 
ically safe to the numerically unsafe region. 

• Eliminate all pentagon tensor integrals by can¬ 
celling terms in the numerator against the propa¬ 
gators wherever possible, after interfering the pen¬ 
tagon amplitude with the Born matrix element. 
The resulting expressions are very large, but nu¬ 
merically stable. 

4. The 0{as) corrections to the Born cross sections 
for the qq and gg initial states due to real gluon emis¬ 
sion, as well as the (g, q)g initiated processes, have been 
computed using both a two cutoff Q and a single cut¬ 
off 1 ^ ^ implementation of the phase space slicing 
(PSS) algorithm. In both PSS methods, the real contri¬ 
bution to the NLO rate is computed analytically below 
the cutoff(s) and numerically above the cutoff(s) and the 
final result is independent of these arbitrary parameters. 
When studying the cutoff dependence of a^Lo, h is cru¬ 
cial to choose the cutoff(s) small enough to justify the 
analytical calculations of the IR divergent contributions 
to d/gtj;, but not so small as to cause numerical instabili¬ 
ties. Finding agreement between the two PSS approaches 
is therefore a strong check of the accuracy of the calcu¬ 
lation. 

In the two cutoff PSS algorithm, the contributions of 
qq, gg tth + g to are first divided into a soft and 
a hard part. 


^ real ^ soft 


' hard 


(5) 


where soft and hard refer to the energy of the final state 
radiated gluon. This division into hard and soft contribu¬ 
tions depends on a soft cutoff, <5^, such that the energy of 
the radiated gluon in the partonic center-of-mass frame is 
considered soft it Eg < ■ The eikonal approximation 

to the soft matrix elements can be taken and the integral 
over the soft degrees of freedom performed analytically. 
Since the (g, q)g initiated process does not develop soft 
singularities, it only contributes to and no soft cut¬ 

off is applied. 

The hard contribution to ij —*■ tth + k (k = g,q,q) 
is further divided into a hard/collinear part, ^]^ard/coiv 
and a hard/non collinear part, aV ,, For gg 

and gg initiated processes the hard/collinear region is 
defined as the region where the energy of the final state 
gluon is Eg > Ss^ and the gluon is radiated from the 
initial massless parton at an angle 9ig , in the ij center-of- 
mass frame, such that (I — cos 0gi) < Sc, for an arbitrary 
small collinear cutoff (5c. The matrix element squared 
in the hard/collinear limit is calculated using the lead¬ 
ing pole approximation, resulting into the convolution of 
the unregulated Altarelli-Parisi splitting functions Pqq or 
Pgg with the corresponding tree level matrix elements 
squared. In the same way, the matrix elements squared 
for the (g, q)g initiated processes are described, in the 
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collinear region (1 — cos0(q_q)i) < 5c, as convolutions of 
the unregulated Altarelli-Parisi splitting functions Pgq 
and Pqg with the qq and gg tree-level matrix elements, 
respectively. The integration over the angular degrees of 
freedom can then be performed analytically. 

The hard gluon emission from the final massive quarks 
never belongs to the hard/collinear region. The contri¬ 
bution from the hard/non collinear region is finite and is 
computed numerically, using standard Monte Carlo inte¬ 
gration techniques. 

Both and o’^^crd depend on the arbitrary cut¬ 

off (5,, and aZrd/coii and hard/non-coU also depend on 
(5c. However, the real hadronic cross section, Ureai, af¬ 
ter mass factorization, is cutoff independent. The cutoff 
independence of the NLO cross section, a^Lo, is shown 
in Fig. 1^, where we assume Sc = ^s/lOO and we let Sg 
vary between 10“^ and 10“^. The hadronic cross sections 
(^soft + c^ hard/coll and <7hard/non-coii include the contribu¬ 
tions of all three tth production channels, qq,gg tthg 
and [q,q)g —> tth{q,q). For (5^ in the range 10“^ —10“^ 
and for 5c = (5^/100, a clear plateau is reached and the 
result is independent of (5^ and Sc- 



Log(8J 

FIG. 1: (7mlo{pp —> tih) calculated in the two cutoff PSS ap¬ 
proach when varying the soft cutoff 5s in the range 
with the collinear cutoff set to (5c = (5s/100, at y^s57=14 TeV, 
for Mh = 120 GeV and g^rrit- The upper plot shows the can¬ 
cellation of the 5s, (5c dependence between a so ft + cr hard/coil 
and <7hard /non-coll- The lower plot shows the dependence of 
onlo on 5s, 5c, with the corresponding statistical errors. 


An alternative method of isolating both soft and 
collinear singularities is to divide the phase space of the 
final state partons into two regions, according to whether 
all partons can be resolved or not. 

The single cutoff PSS technique defines the IR diver¬ 
gent region as that where the final state parton k emitted 
from parton i is not resolved and 

^ik — ‘^Pi * Pk ^ l^min ( 6 ) 


for an arbitrarily small value of the cutoff Smin- The 
partonic real cross section is then written as 


^27 I ^27 -C,Z7 I -27 

= < + ^hard = ^tr +^. 


crossing 


^hard^ 


(7) 



1 10 100 1000 
s.i„ (GeV^) 


FIG. 2: Dependence of aNLo{pp —> tih) on the arbitrary cut¬ 
off Smin of the One cutoff PSS method at 14 TeV, for 

Mh = 120 GeV and g = mt. The upper plot shows the cancel¬ 
lation of the Smin dependence between air = afr + a crossing, 
and ahard- The lower plot shows the dependence of umlo on 
Smin, with the corresponding statistical errors. 


where (t*^ includes both soft and collinear singularities 
and Sr'‘^cird finite. The IR divergent contribution (o'^r^) 
is computed from the crossed process, h —> ijtt+k, where 
ij denote the initial state partons and k = g,q, q. This 
permits a straightforward decomposition of the ampli¬ 
tude in terms of color ordered amplitudes. Using the fac¬ 
torization properties of both the color ordered amplitudes 
and the gluon/quark phase space in the soft/collinear 
limits, the IR divergent contribution can be extracted 
analytically. The factorization of soft and collinear sin¬ 
gularities for color ordered amplitudes has been discussed 
in the literature mainly for the leading color terms, 
0{Nc) l]^, Q. In our case, however, the inclusion 
of the sub-leading terms in 1/Nc is crucial. For the qq 
initial state, these terms were first calculated in Ref. [|j. 

An additional contribution to (T*^ i^^crossing) arises 
when crossing partons i and j back to the initial state, 
due to the mismatch between the collinear gluon radia¬ 
tion from initial and final state partons. distressing ‘con¬ 
tains collinear divergences which are cancelled by the 
parton distribution counterterms when the parton cross 
section is convoluted with the PDF’s. In Fig. |^, the inde¬ 
pendence of CTwlo from the cutoff Smin is demonstrated. 
Again, the hadronic cross sections (t°. -|- Ucrossing and 
(^hard include the contributions from all initial states, 
qq, gg and {q, q)g. Together with Ref. Q, this is the first 
application of the single cutoff PSS approach the calcu¬ 
lation of a cross section involving more than one massive 
particle in the final state. 

The numerical results of both PSS methods agree 
within the statistical errors and within the systematic 
errors of the applied soft and collinear approximations. 
In Ref. 1^, the dipole subtraction formalism has been 
used to extract the IR singularities of the real part. The 
agreement between these three very different treatments 
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of the real IR singularities represents a powerful check of 
the corresponding NLO calculations. 

5. Our numerical results are found using CTEQ4M 
parton distribution functions for the calculation of the 
NLO cross section, and CTEQ4L parton distribution 
functions for the calculation of the lowest order cross sec¬ 
tion [|^. The NLO (LO) cross section is evaluated using 
the 2 (l)-loop evolution of as{fi). The top quark mass is 
taken to be mt = 174 GeV and = 0.116. 

In Fig. ^ we show, for Mh = 120 GeV, the depen¬ 
dence of fJi^o and CTNio on the arbitrary renormaliza¬ 
tion/factorization scale fi = fir = fJ-f- For scales fi larger 
than mt the NLO result is significantly larger than the 
lowest order result. 
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FIG. 3: olo,nlo(jpp tih) as functions of the renormal¬ 
ization/factorization scale fi, at ^/sh = 14 TeV, for Mh = 
120 GeV. 




FIG. 4: a NLO and glo for pp —> tth as functions of Mh, at 
y/sH — 14 TeV, for /r = 2mt -\- Mh . 


Fig. ^ shows both the LO and the NLO total cross 
section for pp tth at = 14 TeV, as functions 
of Mh, for a representative value of the renormaliza¬ 
tion/factorization scale, fj. = 2mt + Mh- We estimate 
the remaining theoretical error, due to the residual fi- 
dependence, to the parton distribution functions, and to 


the experimental error on mt, to be of the order of 15- 
20%. In comparison, the statistical error on the numer¬ 
ical integration is negligible, due to the high statistics 
used in evaluating the total cross section. 


6 . The NLO QGD corrections to the Standard Model 
process pp tth, at = 14 TeV increase the LO 
cross section by a factor of 1.2—1.4 for renormaliza¬ 
tion/factorization scales in the range mt + Mh/2 < fx < 
Amt + 2Mh and Higgs boson masses in the range consid¬ 
ered in this paper. The NLO result shows a drastically 
reduced scale dependence as compared to the Born result 
and leads to increased confidence in predictions based on 
these results. The techniques developed in this calcula¬ 
tion can now be applied to the study of the associated 
bbh production at both the LHC and the Tevatron. 
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